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Abstract 

This work proposes a new model in the context of statistical theory of shape, based 
on the polar decomposition. The non isotropic noncentral elliptical shape distri- 
butions via polar decomposition is derived in the context of zonal polynomials, 
avoiding the invariant polynomials and the open problems for their computation. 
The new polar shape distributions are easily computable and then the inference 
procedure can be studied under exact densities. As an example of the technique, a 
classical application in Biology is studied under three models, the usual Gaussian 
and two non normal Kotz models; the best model is selected by a modified BIC 
criterion, then a test for equality in polar shapes is performed. 
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1 Introduction 



Matrix variate statistical shape analysis has been extensi vely studied in the las t two d ecades 



Dfaz-Garci'a et al. 



Goodalll ( 1991 ) 



Le and Kendall 



(2003)); via affine transformations 



20031 ) , ICaro-Lopera et al I (|2009l ) ) , among many others methods 



(199 



3,, 



by a number of appr oaches: via QR decomposition ( Goodall and Mardia I (119931)); vi a 
SVD decompositio ns 



Diaz-Garcia et al. 



Goodall and Mardia 



(1993) 



(1997), 



Diaz- Garcia et al. 



Drvden and Mardia 



19981) 



and the references there in). However, the polar decomposition has not been included yet 
in the context of shape theory. 

According to the transformation, we say that the shape of an object is all geometrical 
information which remains after filtering out translation, scale, rotation, reflection, uniform 
share, etc, from an original figure comprised in N landmarks in K dimension. Statistical 
shape theory study the mean sha pe of populations in p r esenc e of randomness. 



Some of the classical works ( Goodall and Mardia 



(|1993l )) assume an isotropic Gaus- 



sian model for the landmark matrix in order to obtai n shape densities expanded in known 
polynomials, such as zonal polynomials (|Jamesl (|1964l ). iMuirheadl (|1982l) ); then, generalisa- 
tions for matrix variate sh ape theory under elliptical mo dels appeared via th e SVD method 



( Di'az-Garcia et al. 



([20031 )) and via the affine technique (jCaro-Lopera et al I (|2009l )). How- 
ever, in order to obtain zonal polynomials a partial non isotropy was assumed, otherwise, 
conside r ing a full non isotropy, the densities are expanded in terms of invariant polynomials 



( Davis 



(|l980l )). which are non available for large degrees. 



Now, the isotropic assumption, say & — Ik for an elliptical shape model of the form 

X ~ £nxk{^ x , S x , 0, h), 

restricts substantially the correlations of the landmarks in the figure and it is non appropriate 
for applications. So, we expect the non isotropic model, with any positive definite matrix 
0, as the best model for considering all the possible correlations among the anatomical 
(geometrical o mathematical) points. 

This work solves that problem and sets the non isotropic noncentral elliptical shape dis- 
tributions via polar decomposition in the context of zonal polynomials, avoiding the invariant 
polynomials and the open problems for their computation. The new shape distributions are 
easily computable and then the inference procedure can be studied under exact densities. 

In section [21 the so termed polar shape coordinates are introduced and the main math- 
ematical tools are studied in order to obtained the polar size and shape density. Then the 
polar shape density is derived in section [31 Section 0] studies the central case and the cor- 
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responding invariance under the family of elliptical distributions. Finally, section [5] gives 
explicit densities and performs inference with three models, the Gaussian and two non nor- 
mal Kotz models. 



2 Polar size-and-shape distribution 

Consider a full non isotropy (non singular) elliptical model 

X ~ £ NxK (n x ,5: x ,®,h), 

with generator function h(-). 

In order to avoid the referred problem of invariant polynomials consider the following 
procedure: Let 

X ~ £nxk(p x , S x , 0, h), 
if 1 / 2 is the positive definite sq uare root of the matrix 0, i .e. = (0 1 / 2 ) 2 , with 1 / 2 : 



K x K. 



Gupta and Vargal (|1993l p. 11), and noting that 



X© X X' = X(0- 1/2 0" 1/2 )- 1 X' = X0- 1/2 (X0~ 1/2 )' = ZZ', 



where 



then 



Z = X0~ 1/2 , 



with fi z — /x x -1 ^ 2 , see iGupta and Vargal ( 19931 . p. 20). 

And we arrive at the classical starting point in shape theory where the original landmark 
matrix is replaced by Z = X0 -1 ^ 2 (see Goodall and Mardia (1993), for example). Then 



we can proceed as usual, removing from Z, translation, scale, rotation in order to obtain the 
shape of Z (or X) via QR, SVD, or polar decompositions, for example. 

In this paper we consider a new system of shape coordinates, the polar shape coordinates 
u of X which are constructed as follows: 

LX = Y = RH = rWH = rW(u)H. 

The matrix L is as usual an (N -l)xJV Helmert submatrix and we assume that 
Y~£ N _ lxK (p,i:®I K ,h), ju = L/Lt x , E = LS X L'. 
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Under this approach, Y = RH is the polar decomposition, where R : iV — 1 x Af- lisa 
positive definite matrix and H 6 Vn~i.k- 

It is important to note that only under n = N — 1 the polar approach is valid. 

Recall that for the singular value decomposition Y = PLQ', the polar decomposition of 
Y is given by Y = RH, where R = PL P' and H = P Q'. 

So, we start with a known result, see ICadeti (jl996): 



Lemma 2.1. Let Y : N — 1 x K, then there exist H : N — 1 x N — 1 a positive definite 
matrix and H 6 Vn~i.k such that Y = RH and 

N-l 



{dY) = 11 (Li + Lj)(dR)(HdH'), 



with L = diag(Li, . . . , Ljv-i) an d R = PLP' i.e. L t — Aj(R). 

So the main result of the section follows: 
Theorem 2.1. The polar size- and- shape density is 

N-l 

7T 2 (Li + Lj) 

Proof. The density of Y with f2 = £ V© V is 

/ Y (Y) = — r/i[tr (E^yy + n) ^trE^Y/x'] . 

|S| 2 

If the decomposition Y = RH is performed and the Lemma 12.11 is applied, then the joint 
density of R and H remains 

JV-l 



/ r ,h(R, H) = l<J ,^, K/2 h [tr (S-iR 2 + «) - 2 tr m'S^RH] 



Assum ing that h(-) can be expanded in a convergent power series, see 
il99oL i.e. 



Fang and Zhang 



t=o r ' 



hence 



/k,h(R,H) = ^ /2 ^ ^[tr^R + n)] [tr(-2/j'S _1 RH)]' (HdH'). 
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So, the marginal of R is 



/r(R) 



YliLi+Lj) 

l£ ^hW[tr(E" 1 R+n)] 



E 



/" [tr (HdH') 

** Vw _ i &r 



The integral equals zero when i is odd, then bv I Jamesl (|1964l eq. (22)) 



/ [tr (^m'E^RH)] (HdH') = % £ 7^T C « 



(A) 4* 

Noting that f^ty. ~ tT> so 



JV-l 



/r(R) = 



A v [tr (ir *r 2 + si)] c B (oe^r 2 ) 

2 -iV + i riV _ l[f]|s| f t!( i if ) K 



3 Polar shape density 

Now, observe that R:iV — 1 x N — 1, R>0, contains (JV — l)N/2 different coordinates 
(fy = fjj). Let w(R) the vector consisting of the different elements ry, taken column by 
column. Then the polar shape matrix W, can be written as: 



v(W) = -v(TL), r = ||R|| = VtrR 2 = ||Y|; 
r 



Then bv lMuirheadl ([1982[ ). Theorem 2.1.3, p. 55: 

m m 

(du(W)) = r m J[ sin™- 1 0, /\ d©, A dr, 
which will denoted as 



N(N-l) 

171 = 1, 



i=l 



(dW) - r rn J(u) /\ d0, A dr. 



Thus: 



Theorem 3.1. TTie polar shape density is 

N-l 



2 N-i n ^^ n(A l + A,)J(u) 



/w(W) 



?N-X [f ] |S| 



■EE 



x r 
lo 



(N-l) 2 +2t-l 



/i (2t) [r 2 tr ST *W 2 + tr Si] (dr). 
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Proof. The density of R is 

JV-l 

tv^- H(L i + L 



Let be W(u) = R/r, then the joint density function of W(u) and r is given by 

JV-l 



2 N -^ L ^U(r(X i + X j )) 



/r,w(»)(r,W(u)) = ^ — -g 

» jg> [tr (r'S- 1 W + il)] C K (r^llS- 1 W) 

t=0 K U^Jk 

with to = 7V(7V - l)/2 - 1. Let Xi = A;(W) of W, so if L, = Aj(R), thus L; = rA 4 . Also 
note that: 

1. c K (^ns^w 2 ) = r 2t c K (n-E^w 2 ), 

N-l N-l 

2. IlKA i +A,)=r( JV - 1 )( JV - 2 )/ 2 n(A i + A,), 

3. /^< 2 *) [tr (r 2 S- x W 2 + n)] = fr( 2t ) [r 2 tr E^W 2 + tr O] . 

Collecting powers of r as r m + 2 *+( Ar - 1 )( Ar - 2 )/ 2 = r (JV-i) 2 +2t-i ^ tne ma ,rginal of W is 

2 N-i n &^ TT 1 (A J + A J )J(u) 

f fwl M - c^s-w 2 ) 

X / r (JV-l) a +«-l /l (2t) [ r 2 trS -l W 2 +trn ] (dr) . n 

Jo 

Remark 3.1. Given that H £ Vn-i.k, we cannot classify the polar shape densities by 
including or excluding reflections as in the QR shape distribution cases. 

4 Central case 

The central case of the elliptical polar shape densities follows easily: 
Corollary 4.1. The central polar size- and- shape density is given by 

N-l 

2»-^^R{L l -L,) 

■ WR) = — mew — ft [ ttS ~ lR2 i 
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Proof. Just take fi = in Theorem 12. II and use /i'- -^-) = h(-). 
And finally, we have that 



□ 



Corollary 4.2. The central polar shape density is invariant under the elliptical family and 
it is given by 

2^- V w - 1) ^ w) r[(iv-i) 2 ] 

i<j 



/w(W) 



r*-i [f ] 



JJ(Ai + A,)J(u) (trE- x W 



(N-l)^ 

2\ 5— 



(2) 



Proof. It is straightforward from Theorem 13. II Take (i = 0, and use h^°'(-) = h(-), then 



N-l 



/w(W) 



iW [f ] 1=1 



- J(u) / r^-^-'h [r 2 tr IT *W 2 ] (dr 
Jo 



Let be s = (trE^W 2 ) 2 r, so ds = (tr£ _1 W 2 ) 2 (dr), and 



(JV-1) 2 -1 



(trE^W 2 )- 



if 2\ ^ 



(trE^W 2 ) 2 

(JV-l)" 



trST'W 



JVfJV-l)-!, / 2 



/z (s 2 ) (ds) 



T[(iV-l) 2 ] 



(trE^W 2 ) 2 



Then 



/w(W) = TTTTT^IT II ( A « + J ( u ) (tr= _1 W 



(N-l)' 1 

2\ 2 



2-^+ 2 r JV _ 1 (f ) |s| 



□ 



5 Some particular models 



Finally, we give explicit shapes densities for some elliptical models. 
The Kotz type I model is given by 



Kv) 



n K(N-l)/2 T ( T _ 1 + 



f K(N-l) \ 

' 1 y^expi-Ry}. 



K(N-l) 



So, the corresponding fc-th derivative follows from 
d k 

—^y T ~ 1 exp{-Ry} = 



(-R) k y T - 1 eM-Ry}\l+ £ 



JJ(T-l-i) 



i=0 



(-RyY 



see 



Caro-Lopera e£ 



(2009) 
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It is of interest the normal case, i.e. when T = 1 and R = 5, here the derivation is 
straightforward from the general density. 

The required derivative follows easily, it is, 

K (JV-1) 



and replacing 



/ r (N-lf+2t-l h (2t) [ r 2 trS -l W 2 +tr0 ] 
JO 



dr 



R El^H +2tetr{ Rn} 

f 7V-11 2 



N-l 



2»-^^H(X i + X j )J(u) 



/w(W) 



we have proved that 



EE 



/ r (A r -D 2 +2*-l/,(2*) [ r 2 tr S -1 W 2 + tr ^ 

Jo 



Corollary 5.1. The Gaussian polar shape density is 

N-l 



2 N - 2 J(u) TJ (A; + Xj 



/w(W) 



ctr{-i?fJ} 



r —r N ^ [f ] ( tri?s -i W 2) 



E 



(iV- 



— + *] ^ C K (iJ^S^W 2 ) 



Finally, we propose the result for the Kotz type I model 



fiT _ 1+i£ (iv-i) r / jf(AT-l) 

7 ^ '-^U 1 _i exp{-ity}. 
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Corollary 5.2. The Kotz type I polar shape density is 

JV-l 



2 N - 1 n (k+\j)j(u) 



/w(W) = 



l<3 



EE 



r^_ x [f]|S|f etr{M2}£^ t! 

^ T _ 1+ K(iV-l) +2tr / g(A T_l) ^ 



r T-i 



E^n^-^-^x 



.T-l-i T 



(jV-l) J 



i=0 «=0 



2i?^+^(trIT 1 W 2 ) i ^+ t 



2t 

E 



2t 



m— 1 
oo ^ 2— 1 



m — 1 



n(T i 



E^n^- 1 -") 



(trfi) 



T-l- 



(jV-1) 2 



+ t + i 



i=0 M=0 

Proof. The corresponding fc-th derivative follows from 
^fcJ/ T_1 exp{-ify} = 



2i?^+^(trS- 1 W2)^+' 



(-W^expj-iM 1+E 



n(r-i 



i=0 



(--Ry)" 



(SOO 



Caro-Lopera et 



(2009))), and the corresponding polar shape density is obtained after 
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some simplification like 



N-l 



/w(W) 



t<7 



r*-x [f] 1=1* ^ * ! (K) K 

/•OO 

X / r {N-D 2 +2t-l h {2t) p tr S -1 W 2 + tr „] dr 
•A) 

2 W - X 1] (Ai + A,-)J(u) 



IW [f ] |£|f etr{i?n}^^ 



EE 



TlT-1 

oo i—1 



K(N-l) 
2 



Ein^-^-^x 



(trO 



,T-1- 



(jV-1) 2 



i=0 ti=0 



21 

E 



2/ 



m— 1 
oo ^ i— 1 



ncr-i 



i=0 



E^n^- 1 ---") 



2i? ^ +l+t(trS -i w2) ^)i +t 



(trn) T - x - m -T 



(N-l) 2 



+ t + i 



i=0 u=0 



2J?^+^(trS- 1 W2)^ Mi +* 



□ 



5.1 Example: Mouse Vertebra 



Drvden and Mardia 



(1998) and 



This experiment has been studied in the Gaussian case by 
some references there in. Here we study this data under three different models, the usual 
normal and two non normal Kotz models, the best distribution will determined by applying 
a modified BIC criterion. 

We start with the isotropic version of the Gaussian polar density given in corollary 15. II 

N-l 



J(u) [J (A 4 +A J -)etr(- 



/w(W) = 



MM 

2£ 2 



2 *i^n_iE^p±_ N+2TN _ x r*] ^K(N-l)-(N-ir 

xE- L 



y, C K (^M'WV) 



(-K) 

t=0 K \2 1 ^)k 

Consider now two Kotz models from corollary |5.2l for T = 2 and R = i, we obtain after 
some simplification that 
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N-l 



N Z -K(N-1)-1 



/w(W) = 



J(u) J] (Ai + Xj] 

i<3 



-(N-1) 2 +K(N-1) 



■ etr 



K{N - l)T N -x [f ] V 2S 



E 
E 



(-2t + tr(^S))r[^+t 



(AT-l) 2 



< + 1 



and, for T = 3 and i? = ^ the corresponding non Gaussian isotropic model is given by 



>JV 2_ Jf(w _ 1) + 1 (A ,_ 1) 2_ Jf(jv _ 1) 



N-l 



j( U ) n + \j) 



/w(W) = 

oo 1 

«E5 



*sT(iV - 1)(# (N- l) + 2)rjv_! 



■ etr - 



2£ 2 



t=o 
-4i + 2tr 



2t + 4t 2 - 4< tr 
( fi'fi 



E 



\2S 2 
C K (^/i'W 2 /.) 



\2S 

(N-l) 2 



tr 



V2S 



i + 1 



(iV-1) 2 



+ i + 2 



We contrast these three models via the modified BIC criterion, they will be applied to the 
data of two groups ( small and large) of mouse vertebra, an experiment very detailed in 



Drvden and Mardia 



(119981 ). 

The likelihood based on the exact densities require the computation of the above series, 
a carefully comparison with the known hypcrgcometric of one matrix argument indicates 
that these distribut i ons ca n be obtained by a suitable modification of the algorithms of 
Koev and Edelman I ((20061 ). 

In order to decide which the elliptical model is the best one, different criteria have been 
employed fo r the model selection. We shall consider a modification o f the BIC statis tic as 



discussed in 



Yang and Yangl (|2007l ) , and which was first achieved by iRissanenl (jl978l ) in a 



coding theory framework. The modified BIC is given by: 

BIC* = -2£(/i,S 2 ,/ l ) + n p (log(n + 2) - log 24), 

where £(/x, S , h) is the maximum of the log-likelihood function, n is the sample size and 
Tip is the number of parameters to be estimated for each particular shape density. 
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Table 1: Grades of evidence corresponding to values of the BIC* difference. 



BIC* difference 


Evidence 


0-2 


Weak 


2-6 


Positive 


6-10 


Strong 


> 10 


Very strong 



(1995 


) and 


Rafterv 


(1995) 



criteria have been employed for the model selection. 

In order to apply the above densities we need to restrict the number of landmarks in 
such way that mm(K, N — 1) = N — 1, so in the mouse vertebra data we must select 3 
landmarks of the original 6 points. In the following example we consider the landmarks 1, 
2 and 6 which corresponds to the widest part of the vertebra. 

The maximum likelihood estimators for location and scale parameters associated with 
the small and large groups are summarised in the following table: 

Table 2: Maximum likelihood estimators. 



Group 


BIC* 


Mil 


jJl2 


/U21 


A«22 


a 2 




G 

K : T = 2 














K:T=3 














-129.9719 


-4.5092 


-52.4558 


-18.2941 


3.0990 


32.0686 


Small 


-169.0557 


-6.3553 


-54.5119 


-18.9629 


3.8065 


28.1767 




-139.2333 


-5.6873 


-50.9245 


-20.4403 


4.8584 


27.0289 



Group 


BIC* 


Mil 


Ml2 


A«21 


^22 


a 2 




G 

K : T = 3 














K :T=2 














-128.2599 


-28.0776 


-73.1440 


-32.0698 


13.1131 


75.7233 


Large 


-204.2375 


-20.1702 


-95.3134 


-41.9462 


9.8279 


84.3441 




-167.2111 


-26.0985 


-82.8228 


-36.3675 


12.3299 


75.0727 



According to the modified BIC criterion, the Kotz model with parameters T = 2, R = \ 
and s = 1 is the most appropriate for the small group, instead the Kotz distribution with 
parameters T = 3, R = j an d s = 1 models the large group. There is a very strong difference 
between these models and the classical normal in this experiment. 

Let fi 1 and fj, 2 be the mean polar shape of the small and large groups, respectively. We 
test equal mean shape under the best models, and the likelihood ratio (based on —2 log A « 
xl) for the test Hq : ii x = fj, 2 vs H a : fij ^ /x 2 , provides the p- value 0.99, which means that 
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there extremely evidence that the mean shapes of the two groups are equal. 

For any elliptical model we can obtain the polar shape density, however a nontrivial 
problem appears, the 2i-th derivative of the generator model, which can be seen as a partition 
theory problem. For The general case of a Kotz model (s ^ 1), and another models like 
Pearson II and VII, Besse l, Jensen-logistic, we can use formulae for these derivatives given by 
Caro-Lopera et al I (|2009f) . The resulting densities have again a form of a generalised series 



of zonal polynomials which can be computed efficiently after s ome modification of existing 



works for hypergeometric series (see 



Koev and Edelman 



(2006)), thus the inference over an 



exact density can be performed, avoiding the use of any asympt o tic di stribution, and the 



initial transformation avoids the invariant polynomials of 
seems can not be computable for large degrees. 



Davis 



(1980), which at present 
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